Cosmological constraints in the presence of ionizing and resonance radiation at 

recombination 
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With the recent measurement of full sky cosmic microwave background polarization from WMAP, 
key cosmological degeneracies have been broken, allowing tighter constraints to be placed on cos- 
mological parameters inferred assuming a standard recombination scenario. Here we consider the 
effect on cosmological constraints if additional ionizing and resonance radiation sources are present 
at recombination. We find that the new CMB data significantly improve the constraints on the 
additional radiation sources, with log 10 [e a ] < —0.5 and log 10 [e;] < —2.4 at 95% c.l. for resonance 
and ionizing sources respectively. Including the generalized recombination scenario, however, we 
find that the constraints on the scalar spectral index n 3 are weakened to n 3 — 0.98 ± 0.03, with 
the n B = 1 case now well inside the 95% c.l.. The relaxation of constraints on tensor modes, scale 
invariance, dark energy and neutrino masses are also discussed. 



I. INTRODUCTION 

The recent measurements of the Cosmic Microwave 
Background (CMB) flux provided by the three year 
Wilkinson Microwave Anisotropy Probe (WMAP) 
mission (see 0, 0, i || have confirmed several of the 
results already presented in the earlier data release, 
but also pointed towards new conclusions. The bet- 
ter treatment of systematics in large scale polarization 
data, in particular, has now provided a lower value for 
the optical depth parameter r. This, together with 
an improved signal in the temperature data at higher 
multipoles, has resulted in a lower value of the spec- 
tral index parameter n s — 0.959 ± 0.016. A deter- 
mination of this parameter can play a crucial role in 
the study of inflation. Soon after the WMAP data 
release, several papers have indeed investigated the 
possibility of discriminating between single-field infla- 
tionary models by making use of this new, high qual- 
ity, dataset 0, i, 0, S, S, 03, EE El El. One of 
the main conclusions of these papers is that some in- 
flationary models, such as quartic chaotic models of 
the form V((j>) ~ A<fi 4 , may be considered ruled out 
by the current data while others, such as chaotic in- 
flation with a quadratic potential V(<j>) ~ m 2 (j> 2 are 
consistent with all data sets. 

While the WMAP result is of great importance 
for inflationary model building, one should be care- 
ful in taking any conclusion as definitive since the 
constraints on n s are obtained in an indirect way 
and are, therefore, model dependent. Similar con- 
siderations applies to other cosmological constraints, 
such those on the dark energy equation of state and 
neutrino masses. Combining CMB anisotropies with 
galaxy clustering and supernovae type la data, the 
dark energy equation of state parameter (dark en- 
ergy pressure over density) has been constrained to 
w = -1.08 ± 0.12 at 95% c.l. (see Q). Using the 
same dataset, but under the assumption of a cosmo- 
logical constant, it is possible to constrain the neu- 



trino masses to 53 m * < 0.66ey at 95% c.l. where 
i = 1,..,3 and indicates the neutrino flavor. Again, 
while those constraints play a very important role in 
our understanding of the dark energy component and 
neutrino physics, they are obtained in an indirect way 
and under several assumptions. 

The importance of the model dependency of the cos- 
mological constraints has been recently discussed by 
several authors. The impact of isocurvature modes 
on the determination of the neutrino mass 1 141 . dark 
energy properti es 11511 . scalar spectral index [lj| and 
baryon density [17| is just one example. 

Here we investigate possible deviations in the mech- 
anism on which CMB anisotropies are highly depen- 
dent: the process of recombination. 

In a previous paper [l8j . we analyzed modified re- 
combination processes in light of the WMAP first year 
data. Here we assess the improvements given by more 
recent data, in particular the inclusion of CMB po- 
larization spectra, and also extend the analysis to a 
larger set of parameters. We will indeed not only pro- 
vide new and more stringent constraints on modified 
recombination but we also consider its impact on in- 
flationary, dark energy and neutrino parameters, 

The recombination process can be modified in sev- 
eral ways. For example, one could use a model- 
independent, phenomenological approach such as in 
[l9| where models are specified by the position and 
width of the recombination surface in redshift space. 
Here we instead focus on theoretically motivated 
mechanisms based on extra sources of ionizing and 
resonance radiation at recombination (see e.g. (20|). 
While the method we adopt will be general enough 
to cover most of the models of this kind, as discussed 
in the next section, we remind the reader that there 
exist other ways in which to modify recombination, 
for instance, by having a time- varying fine-structure 
constant fp^O]!. 

Following the seminal papers [13, detailing 
the recombination process, further refinements to the 
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standard scenario were developed [24|, allowing pre- 
dictions at the accuracy level found in data from 
the WMAP satellite and the future Planck satellite 
[25l [26| . With this level of accuracy, it becomes 
conceivable that deviations from standard recombina- 
tion maybe be detectable [2(| H3, HH , although fur- 
ther refinements could be required to get the Thom- 
son visibility function below percent level accuracy 

sunn. 

The paper proceeds as follows: in section [IT] we de- 
scribe a model which can produce deviations from the 
standard recombination scenario. In IIII1 we describe 
how these deviations might affect the CMB temper- 
ature and polarization power spectra and conduct a 
likelihood analysis using the recent CMB data from 
WMAP and other cosmological observables. In par- 
ticular, we will study the impact that a modified re- 
combination scheme can have on several cosmological 
and astrophysical parameters. In llVl we draw together 
the implications of the analysis. 



II. A MODIFIED IONIZATION HISTORY 

The evolution of the ionization fraction, x e , of 
atoms, number density n, can be modeled in a sim- 
plified manner for the recombination of hydrogen, 
Ml, 



dx e 
~~dt 



std 



c 



a c nx\ — b c (l — x e ) exp 



k B Tj 



(1) 



where a c and b c are the effective recombination and 
photo-ionization rates for principle quantum numbers 
> 2, AB is the difference in binding energy between 
the I s * and 2 nd energy levels and 



C = 



1 + KA ls2s n ls 
1 + K(A ls2s + b c )n u 



K = 



Ail 



8nH(z) 



(2) 



where A Q is the wavelength of the single Ly-a transi- 
tion from the 2p level, Ai S 2 S is the decay rate of the 
metastable 2s level, ni s = n(l — x e ) is the number of 
neutral ground state H atoms, and H (z) is the Hubble 
expansion factor at a redshift z. 

We include the possibility of extra photons at key 
wavelengths that would modify this recombination 
picture, namely, resonance (Ly-a) photons with num- 
ber density, n a , which promote electrons to the 2p 
level, and ionizing photons, n, ■., [ill. I20L I27L I2H] 



— = e a (z)H(z)n, — 



(z)H(z)n. (3) 



which leads to a modified evolution of the ionization 
fraction 



dx 



-£ = - jt L« - C ^ H ~ (1 - °)^ H - ( 4 ) 



dt 



Extra photon sources can be generated by a variety 
of mechanisms. A widely considered process is the 



decay or annihilation of massive particles [2fJ, |32|, |33j, 
H3, [H, [H, H3, HH . The decay channel depends on the 
nature of the particles, and could, for example, include 
charged and neutral leptons, quarks or gauge bosons. 
These particles may then decay further, leading to a 
shower /cascade that could, amongst other products, 
generate a bath of lower energy photons that could 
interact with the primordial gas and cosmic microwave 
background. Interestingly these models, as well as 
injecting energy at recombination, z ~ 1000, boost the 
ionization fraction after recombination and can distort 
the ionization history of the universe at even later 
times, du ring galaxy formation and reionization z ~ 
5 - 10 d HI SO, El S3. 0ther mechanisms include 
evaporation of black holes [U S3 or inhomogenities 
in baryonic matter (27j. 

We employ the widely used RECFAST code [Hj|, in 
the cosmomc package [45| modifying the code as in (f?|) 
to include two extra constant parameters, e a and q. 
In addition to the ionizing sources, we assume a single, 
swift reionization epoch at a redshift z re . 

In Figure [T] we show the effect of additional res- 
onance and ionizing radiation on the CMB TT, TE 
and EE spectra, in comparison to a fiducial best fit 
model to the WMAP 3-year data. From identical ini- 
tial power spectra, the inclusion of additional reso- 
nance photons slighly boosts the ionization fraction 
at and after recombination, suppressing TT power at 
small scales, while the large scale EE spectra is largely 
unaffected. Ionizing photons significantly boost the 
ionization fraction post recombination and therefore 
as well as significantly suppressing TT power on small 
scales, they also can generate a boost in the large scale 
EE signal akin to an early partial reionization. 



III. LIKELIHOOD ANALYSIS 

The method we adopt is based on the publicly avail- 
able Markov Chain Monte Carlo package cosmomc 
?45|. We sample the following dimensional set of cos- 
mological parameters, adopting flat priors on them: 
the physical baryon and Cold Dark Matter (CDM) 
densities, = flbh 2 and w c = fl c h 2 , the ratio of the 
sound horizon to the angular diameter distance at de- 
coupling, 9 S , the scalar spectral index, n s , and the 
optical depth to reionization, r. As described in the 
previous section, we modify recombination by consid- 
ering variations in the e a and £j parameters. Further- 
more, we consider purely adiabatic initial conditions 
and we impose flatness. We also consider the possi- 
bility of having a tensor (gravity waves) component 
with amplitude r respect to scalar, a running of the 
spectral index dn s /dlnk at k = {).{){)2h~ 1 Mpc and a 
non-zero, degenerate, neutrino mass of energy density: 



n v h 2 = 



92.5ey 



(5) 
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FIG. 1: (From left to right) The comparison of ionization fraction evolution, and TT (left) ,TE (center) and EE (right) 
CMB spectra comparing a best fit ACDM fiducial model to models with the same cosmological parameters but with 
additional resonance (top) and ionizing (bottom) radiation. WMAP binned data are shown as blue shaded regions. 



Finally, we will also investigate the possibility of a 
dark energy equation of state, w, different from —1 
but constant with redshift. The MCMC convergence 
diagnostics is done on 7 chains though the Gelman 
and Rubin "variance of chain mean" / "mean of chain 
variances" R statistic for each parameter. Our ID 
and 2D constraints are obtained after marginaliza- 
tion over the remaining "nuisance" parameters, again 
using the programs included in the cosmomc pack- 
age. In addition to the WMAP data, we also consider 
the constraints on the real-space power spectrum of 
galaxies from the Sloan Digital Sky Survey (SDSS) 
[4q | . We restrict the analysis to a range of scales over 
which the fluctuations are assumed to be in the linear 
regime (k < 0.2/i _1 Mpc). When combining the mat- 
ter power spectrum with CMB data, we marginalize 
over a bias b considered as an additional nuisance pa- 
rameter. Furthermore, we make use of the Hubble 
Space Telescope (HST) measurement of the Hubble 
parameter Hq = 100ft km s _1 Mpc _1 [13] by multi- 
plying the likelihood by a Gaussian likelihood function 
centered around h = 0.72 and with a standard devia- 
tion a = 0.08. When considering dark energy models, 
we also include information from luminosity distance 
measurements of type la Supernovae from the recent 
analysis of [48| . Finally, we include a top- hat prior on 
the age of the universe: 10 < to < 20 Gyrs. 



IV. RESULTS 

Our main results are plotted in Figure [2] where we 
show the 68% and 95% c.l. on the n s — log 10 (e Q ), 
<7 8 -log 10 (e a ), n s -log 10 (ei) and cr 8 - log 10 (e;) plane. 



In the top portion of Figure [2j we consider only the 
WMAP data (plus a prior on the Hubble parameter), 
while in the lower portion, we add SDSS. Let us first 
consider the case of WMAP alone. As we see, using 
this dataset alone, we can put interesting new bounds 
on the recombination parameters. Marginalizing over 
the remaining, "nuisance" , parameters we indeed ob- 
tain log 10 (e a ) < -0.81 and log 10 (ej) < -2.31 at 95% 
c.l.. 

As suggested by Figure [TJ we find ionizing photons 
are better constrained with current data since the ion- 
ization fraction is significantly boosted at and beyond 
the onset of recombination. This results in a suppres- 
sion of TT power and boosting of EE power even on 
large scales, well constrained by WMAP data. Reso- 
nance photons have a more subtle effect only slightly 
increasing the ionization fraction after the onset of 
recombination. This leads to a suppression of small 
scale TT power but little effect on large scale EE. 
The constraints on both types of radiation should be 
noticeably improved therefore by future experiments, 
such as the planned PLANCK satellite, which better 
characterize small scale TT and EE power [28| . 

Moreover, there is a clear degeneracy between 
log 10 (e Q ) and the spectral index n s . Indeed, a modifi- 
cation of the recombination scheme allows us to sup- 
press the amplitude of the peaks in the CMB power 
spectrum in a way similar to a later recombination but 
without altering the large-scale polarization signal. 
This changes in a drastic way the constraints on the 
scalar spectral index and og. Marginalizing over the 
recombination parameters, we get n s — 0.978^Q Q2g 
and a 8 = 0.80±g;g| at 95% c.l.. Those results should 
be compared with the constraints n s — 0.9591q'q27 
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FIG. 2: The 68% and 95% likelihood contours in the n s 
and (J8 vs e a plane (left) and n s and as vs Ej (right). 
The analysis considers (top/blue) the 3-years WMAP data 
and a HST prior on the Hubble parameter, h, alone and 
(bottom/red) also including SDSS galaxy matter power 
spectrum data. 
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FIG. 3: The effect of a modified recombination scheme 
on inflationary parameters in a WMAP+SDSS analysis. 
In the top panel we plot the constraints in the n s — r 
plane. The filled contours (68% and 95%) are obtained af- 
ter marginalization over extended recombination parame- 
ters while the empty contours assume standard recombina- 
tion. In the bottom panel, we show the 68% and 95% likeli- 
hood contours in the dn„/dlnk vs e a plane. Non-standard 
recombination shifts the ~ ler evidence for running in the 
standard case to a null result. Future evidence for run- 
ning may be interpreted as evidence for a non-standard 
recombination scheme. 



and erg = 0.781q q7, again at 95% c.L, obtained using 
the same dataset and priors but with standard recom- 
bination. 

Including SDSS data, as is shown in the lower panel 
in Figure [2j does not significantly improve our con- 
straints on e a and e^. The SDSS power spectrum in- 
deed prefers a higher value of the as parameter than 
WMAP. While the tension is not strong enough to 
provide any evidence for modified recombination, the 
constraints are lowered to log 10 (e Q ) < —0.51 for e a 
and almost stable to log 10 (ei) < —2.24, for at 
95% c.L. The constraints on n s and as are also af- 
fected. Including SDSS we find n s = 0.994±£^ and 
as = 0.87to° o l at 95% c.L. 

Interestingly, we find that the constraints on other 
key parameters (r or fij,) are robust to the modifica- 
tions in the recombination scenario. It is interesting 
to extend the analysis to other inflationary parame- 



ters such as the amplitude of a tensor component r 
or a running of the spectra index dn s /dlnk. In Fig- 
ure[3](top panel) we plot the 68% and 95% likelihood 
contours in the n s — r plane in the standard and in 
the generalized recombination case. As one can see, 
relaxing our knowledge about recombination strongly 
affects the final constraints: the scalar spectral index 
can be more consistent with n s > 1 and the upper 
limit on the tensor component can be a factor 2 larger 
than in the standard case. As one can see from the 
bottom panel of Figure[3]a degeneracy between e a and 
the running dn s /dlnk is also present. Standard anal- 
yses prefer a negative running of the spectral index 
with significance slightly above la (see [l|). This can 
be compensated for by a non-standard recombination 
with e a > 0.1. 

In Figure [H we report on the impact of non- 
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bined, indication for w < — 1 could, therefore, provide 
a hint of a non-standard recombination process and 
one should be careful in interpreting it as evidence for 
a phantom- like dark energy component. In a more 
generalized recombination scenario, we find the con- 
straints on w are relaxed to w = — 1.241°' 44 at 95% 
c.L. 

Finally, in Figure [5J we report the constraints on 
neutrino masses. As one can see, non-standard re- 
combination also relaxes constraints on this parame- 
ter. We find that values as large as Sm^ ~ 1.2eV 
are consistent with the data, relaxing by ~ 50% the 
standard constraint £m„ < 0.72eV (see e.g. [H, E§|). 



FIG. 4: The impact of a modified recombination scheme 
on constraining a constant dark energy equation of state, 
w. We show the 68% and 95% likelihood contours in 
the w vs e a plane from WMAP+SDSS+HST+SN-la (see 
text). Non-standard recombination relaxes the constraints 
towards more negative values for w 
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FIG. 5: The effect of a modified recombination scheme 
on constraining neutrino masses. We show the 68% and 
95% likelihood contours in the Em„ vs e a plane from 
WMAP+SDSS+HST (see text). Non-standard recombi- 
nation relax the constraints towards larger masses 



standard recombination on the equation of state pa- 
rameter, w. We find an important degeneracy only 
with e Q ; allowing e a to vary enlarges the constraints 
on w towards more negative values. A future, com- 



V. CONCLUSIONS 

In this paper, we update the upper bounds that 
can be placed on the contribution of extra Ly-a 
and ionizing photon-producing sources in light of the 
new WMAP data. We find that, adopting a simple 
parametrization using constant effective values for e a 
and Si, the WMAP data constraints log 10 [e Q ] < —0.5 
and log 10 [£i] < —2.4 at the 95% level. Physically mo- 
tivated models for non-standard recombination which 
generate ionizing and resonance radiation, like those 
based on primordial black hole or super-heavy dark 
matter decay, remain feasible. 

We find that a modified recombination scheme may 
affect the current WMAP constraints on inflationary 
parameters like the spectral index n s and its running. 
In particular, if recombination is changed, Harrison- 
Zel'dovich spectra with n s = 1, larger tensor modes 
and positive running are in agreement with observa- 
tions. Moreover, constraints on particle physics pa- 
rameters like the neutrino mass are also relaxed when 
non-standard recombination is considered. 

Future observations in both temperature and po- 
larization, such as those expected from the Planck 
satellite [28[ , will provide more precise small scale TT 
and EE measurements needed to more stringently test 
these models and, crucially, will reduce the depen- 
dency of other cosmological parameters on them. 
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